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Abstract 

Recent observations suggest that blobs of matter are ejected with ultra- 
relativistic speeds in various astrophysical phenomena such as supernova ex- 
plosions, quasars, and microquasars. In this paper we analyze the gravita- 
tional radiation emitted when such an ultra-relativistic blob is ejected from a 
massive object. We express the gravitational wave by the metric perturbation 
in the transverse-traceless gauge, and calculate its amplitude and angular de- 
pendence. We find that in the ultra-relativistic limit the gravitational wave 
has a wide angular distribution, like 1 + cos#. The typical burst's frequency 
is Doppler shifted, with the blue-shift factor being strongly beamed in the 
forward direction. As a consequence, the energy flux carried by the gravita- 
tional radiation is beamed. In the second part of the paper we estimate the 
anticipated detection rate of such bursts by a gravitational-wave detector, for 
blobs ejected in supernova explosions. Dar and De Rujula recently proposed 
that ultra-relativistic blobs ejected from the central core in supernova explo- 
sions constitute the source of Gamma-ray bursts. Substituting the most likely 
values of the parameters as suggested by their model, we obtain an estimated 
detection rate of about 1 per year by the advanced LIGO-II detector. 



I. INTRODUCTION 



Relativistic jets seem to be emitted by astrophysical systems wherein mass is accreted 
at high rate from a disk to a central compact object (for a review see []IJ). Astrophysical 
observations suggest that blobs of plasma are ejected with ultra-relativistic velocities in 
supernova explosions @, microquasars GRS1915+105 §, §, [§, § and GRO J165-40 [§, 
and in active galactic nuclei hosting a massive black hole. 

Recently Dar and De Rujula proposed a new model for the origin of Gamma ray bursts, 
the cannonball model H, in which the bursts are sourced by ultra-relativistic blobs of matter 
emitted in supernova explosions. According to this model, in a typical supernova explosion 
some portion of the expanding mass falls back on the central core and forms an accretion 
disc. Accreted matter is then ejected as blobs of plasma in the polar directions, with a 
large Lorenz factor 7 of order 10 3 . A strongly beamed burst of Gamma ray is created 
when such a blob hits a shell of expanding matter ejected earlier in the supernova explosion 
process. According to the cannonball model, in a typical supernova explosion a few such 
ultra-relativistic blobs are ejected in each of the two polar directions, yielding a Gamma ray 
burst made of a few strong peaks. 

In this paper we investigate the gravitational radiation from such ultra-relativistic blobs 
of matter, and evaluate the anticipated rate of detections by the gravitational-wave detector 
LIGO. Since the motion is relativistic the quadrupole formula cannot be used in this prob- 
lem. Instead, we solve the linearized Einstein equations using the Lienard-Wiechert formula 
(generalized to the gravitational case). 

Gravitational radiation is emitted whenever the blob changes its velocity (the gravita- 
tional field involved in a motion with constant velocity is non-radiative). In this paper we 
focus on the radiation emitted when the blob is ejected from the central object and is ac- 
celerated to a large Lorenz factor. A burst of gravitational radiation is also emitted when 
the blob hits the ejecta, but this burst appears to be weaker by several orders of magnitude 
and we shall not consider it here. 
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The analysis throughout this paper is essentially free of assumptions about the values 
of the astrophysical parameters involved. We merely assume that (i) a blob of matter is 
ejected from a massive object (a "star") and is accelerated to a Lorenz factor 7 >> 1, 
and (ii) the blob's energy 7m is small compared to the star's mass M. (Both assumptions 
7 >> 1 and 7m « M are not necessary for the analysis, but they significantly simplify it.) 
We obtain a general expression for the amplitude of the gravitational wave as a function of 
direction. We also calculate the directional dependence of the gravitational wave's observed 
frequency. These expressions involve two astrophysical parameters: the blob's energy 7m 
(for the amplitude), and the typical time scale At of acceleration (for the frequency). We 
then derive a general expression for the anticipated detection rate. This expression depends 
on three more parameters: the event rate per unit volume (e.g. the event of supernova 
explosions), the detector's sensitivity, and the detector's optimal frequency. Substituting the 
most likely values of the astrophysical parameters, as suggested by the cannonball model, we 
obtain an anticipated detection rate of about 1 per year by the advanced LIGO-II detector. 
This detection rate is not certain, however, because of an uncertainty in the astrophysical 
parameters involved. In particular, the detection rate is proportional to a third power of 
the blob's energy, and the uncertainty in the latter may change the detection rate by one or 
two orders of magnitude. 

Our analysis shows that despite the large Lorenz factor, the gravitational field, expressed 
in terms of the metric perturbation, is not strongly beamed in the forward direction. Rather, 
at the ultra-relativistic limit the directional dependence of the transverse-traceless (TT) 
metric perturbation is like 1 + cos 9, where 9 is the angle between the particle's velocity 
and the spatial direction vector from the source point to the observer. Thus, whereas the 
Gamma-ray burst can only be observed in a very small solid angle comparable to 7~ 2 (due 
to the strong beaming of the electromagnetic radiation), the gravitational signal may be 
observed in a wide solid angle, effectively ~ 2ir steradians (unless At is too large - see 
section VI). On the other hand, the observed frequency is strongly Doppler blue-shifted in 
the forward direction. As a consequence, the energy flux carried by the gravitational waves 
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is beamed in the forward direction. 

The gravitational radiation emitted in this process has two special features which dis- 
tinguish it from most other sources. First, we are dealing here with a "burst with memory" 
||. That is, at the end of the process (e.g. ejection of a blob) the metric perturbation 
amplitude does not return to its original value (see section IV). Secondly, according to the 
cannonball model, in a typical supernova event several blobs are emitted (in each of the two 
polar directions). Consequently, the gravitational signal will be composed of a few separate 
bursts. Therefore, once such an event is detected, it may be easy to distinguish it from other 
sources of gravitational radiation. 

In section II we calculate the metric perturbation produced by the moving blob. We 
first calculate it in the Lorenz gauge, using the gravitational analog of the Lienard-Wiechert 
formula. Then we transform the metric perturbation to the TT gauge. In section III we 
obtain the angular dependence of the wave's amplitude in the ultra-relativistic limit. We 
find this amplitude to be proportional to 1 + cos# (with a narrow "hole" at the center, 
whose angular width is ~ We also discuss the relation between our results and a 

previous work by Dray and tHooft (DtH) ||, who analyzed the gravitational field of a 
particle moving at a (fixed) ultra-relativistic speed. Then in section IV we calculate the 
total change in the metric perturbation that occurs when a massive star emits an ultra- 
relativistic blob. We show that this change is non- vanishing (namely, this is a "burst with 
memory"). Furthermore, this change is (at the leading order) equal to the contribution of 
the blob itself to the metric perturbation. In section V we obtain the angular dependence 
of the observed burst's frequency, which is strongly Doppler blue-shifted in the forward 
direction. We show that unlike the metric perturbation, the energy flux is indeed beamed 
in the forward direction, 9 ~ I/7. 

In section VI we derive the expression for the anticipated detection rate, as a function 
of the various parameters involved. Finally, in section VII we substitute in this general 
expression the astrophysical parameters emerging from the cannonball model, as well as the 
LIGO-II detector's parameters. Of these parameters, two have the largest uncertainty: the 
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blob's energy E = ^mc 2 and the acceleration time At. The cannonball model yields an 
order-of-magnitude estimate for the blob's energy: E ~ 10 52 erg. The parameter At has a 
larger uncertainty; however, this parameter does not affect the detection rate as long as it 
is smaller than the detector's typical time scale. And, even if it is larger, it only affects the 
detection rate through its first (inverse) power. Therefore, the main uncertainty seems to 
come from the energy parameter, which enters the detection rate as E 3 . With the above 
value for E, we obtain a detection rate of about 1 event per year in the advanced LIGO-II 
detector (provided that At is not too large). 

For the above value of E, the maximal distance for observation by LIGO-II is found to 
be -R max ~ 15Mpc. Since this corresponds to z « 1, we ignore the cosmological redshift 
effects throughout the paper. 

Several authors, mostly during the 1970's, investigated the gravitational radiation emit- 
ted from ultra-relativistic sources ( [jlOfl , [pTT[] , [fL2] , fl3fl , [ |14|j , ||15|| , |T6|| and references therein). 



These authors considered a variety of model problems, used several methods of calculation, 
and studied various aspects of the gravitational field emitted. We haven't found any previ- 
ous work which covers the problem that concerns us here. The closest we have found is the 
analysis by Adler and Zeks who considered a similar problem of a supernova explosion. 
However, they calculated the TT wave-form only in the special case of two equal masses, 
whereas the astrophysical situation that concerns us here is m « M. (Indeed, in the 
special case of two equal-mass blobs ejected simultaneously in the two polar directions, our 
result agrees with Ref. |12[ - see section IV.) Peters [|1(J, and later Kovacs and Thorne |13 



15], considered ultra-relativistic encounters, but their analysis is restricted to the case of a 



large impact parameter, whereas in our problem the impact parameter vanishes. Ruffini [5 I 



Smarr |TJj and D'eath |]T6| studied the ultra-relativistic head-on collision of two black holes 



(or of a small object with a massive black hole |TT||). This situation differs from our case, in 
which the star and the blob are both weak-field objects. (There also is a difference in the 
time direction, i.e. an explosion instead of a collision, but at least in our star-blob system 
the gravitational radiation does not care about this change in the direction of time - see 



section V). Naively one might expect that the gravitational waves produced in the collision 
will not be sensitive to the nature of the objects involved, as long as the latters are small. 
Our analysis, however, suggest the contrary for the head-on case. In our star-blob problem 
we find the energy flux of gravitational radiation to be beamed in the forward direction (see 
section V). No such beaming occurs in the head-on collision of two black holes [13 



This difference between the two problems has a simple intuitive reason, which we discuss in 
section V. 

We use the signature ( h ++)• Since most of the paper deals with basic general- 

relativistic analysis, we use general-relativistic units G = C = 1. Only in section VII, in 
which we put astrophysical numbers, we retain the values of G and C in standard physical 
units. 

II. GRAVITATIONAL FIELD OF A RELATIVISTIC PARTICLE 

For a weak gravitational field, the linearized Einstein equations (expressed in Cartesian 



coordinates x a ) read fl? 



IGttT = -h a -n h a al3 +h a +h a (1) 

where h^ v = h^ v — \r]^ v h\ and h^ u is the metric perturbation. Under the Lorenz gauge 
conditions h m a = 0, equation ([!]) reduces to 

V^=-167rT^. (2) 

Consider a point mass m (a "particle") moving along a world line r a (r), where r is the 
proper time and r a denotes the particle's location in Cartesian coordinates. The energy 
momentum tensor of such a point mass is given by 

T afi (x)=m [ u a (r) v? (r) 5 (4) [x - r (r)] dr , (3) 



where u a = dr a /dr is the particle's 4- velocity. The retarded solution of Equation (|2]) for 
such a source term is obtained by a straightforward generalization of the Lienard-Wiechert 
potentials: 
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-pa/3 , , , m q (r) it^ (r) 
h (x) = 4m ^ 



(4) 

T=To 



-m 7 • [x — r (r)] 

This expression is to be evaluated at the retarded time tq, which is the intersection time of 
r a (r) and the observer's past light-cone. The metric perturbation in the Lorenz gauge is 
then given by 



u a (r) vP (r) + K] al3 



(5) 



2'' ' 7 -u 7 • [x - r (r)] 7 

In the next step we transform from the Lorenz gauge to the TT gauge, which is best 
suited for calculating the response of a gravitational- wave detector. The metric perturbation 
in the TT gauge, which we denote hJJ, includes only space-space components, namely 
hfj" = 0. This spatial part is obtained from h^ u by fl7| 



h TT = PhP --P -Tr(hP) . (6) 
2 

Here h TT , h, and P are 3x3 spatial matrices, where h represents the spatial part of h^, 
and P is a projection operator defined by Pjk = Sjk — fiju^, where n is the unit spatial 
direction vector from the (retarded) source point to the observer. Hereafter Latin indices 
run over the three spatial components. [In Eq. (Q) we have omitted the indices for brevity, 
and we use the standard matrix product notation]. 

Before proceeding with the calculation, there is a subtlety that must be addressed. The 
projection operation (^) applied to a metric perturbation h constitutes a gauge transforma- 
tion only if h is a pure gravitational-radiation field. In our case, for a particle moving at a 
fixed speed, h is non-radiative. Despite of this, the application of Eq. ) to our problem is 
justified, because of the following reason: The physical quantity that will concern us in this 
paper is not the value of h, but rather the change in h that occurs during an astrophysical 
process. This change, which we denote Ah, occurs when the particle changes its velocity 
due to interaction with another object. The quantity Ah represents a pure radiation field, 
and hence applying the projection (§) to it yields a valid gauge transformation. It is this 
radiative piece Ah (in fact, its TT part) that is relevant for detection over astrophysical 



distances. Obviously, in order to have a nonvanishing Ah we must consider a system of two 
particles (or more), interacting with each other. Then we have to sum over the contributions 
to h from all components of the system. Correspondingly, the quantity relevant for gravi- 
tational wave detection is the change in this sum, namely Ah = A(Eh), where £ denotes a 
summation over the components of the system. This quantity is a pure radiation field, and 
from the linearity of the projection @, the TT-part of Ah is given by 

Ah TT = (Ah) TT = A(£(/i TT )) . (7) 

In what follows we shall calculate h TT (for a single object), and in section IV we shall 
construct from it the quantity Ah TT for the situation of a blob ejected from a star. (In fact, 
we shall show that at the leading order Ah TT is nothing but h TT of the ejected blob.) 

Proceeding with the calculation of h TT , one finds that the term proportional to r] a ^ in 
Eq. (H), which represents a pure trace term, vanishes upon the projection (|6]), therefore 

h TT = PhP -ip -Tr(hP) , (8) 

where h denotes the 3x3 spatial part of h^ v . It is useful to decompose this expression into an 
amplitude factor and a directional factor. We define the particle's 3- velocity, v = df/dt, the 
particle's speed v = \v\, and the unit 3-vector in the velocity direction, v = v/v (hereafter a 
bar denotes a spatial 3-vector, and a unit 3-vector is denoted by a hat). Defining Wij = ViVj, 
we find hij = h Wij where 

= 4m(/? 7 ) 2 = 4 7 m/3 2 
-m 7 • [x - r (r)] 7 i2(l-/?cos0) ' 1 ' 

Here j3 and 7 denote the standard special-relativistic quantities (3 = v /c and 7 = (1— /3 2 ) -1 / 2 , 

R is the spatial distance between the evaluation point x a and the source point r a , and 9 is 

the angle between n and v, i.e. n ■ v = cos 9. Equation (p|) now reads 



h TT = he 



PwP -^P-Tr {wP) 



h r 



PwP - -Psin 2 9 
2 



(10) 



This expression enfolds the information on both the amplitude and the polarization of the 
emitted gravitational wave. The direction of polarization is dictated by the transversality 
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condition and the direction of motion. Thus, at a given space-time point, the maximal 
detection amplitude, which we denote h + , is achieved for a detector's arm directed perpen- 
dicular to n in the h-v plane. The same amplitude (but with opposite sign) is obtained 
in the perpendicular transverse direction, i.e. in the direction perpendicular to both n and 
v, which we denote n + . The gravitational wave's amplitude h + can thus be obtained by 
projecting h TT on the direction n + : 

h + = -nfhJ?n+. (11) 

(Hereafter, a repeated spatial sub-index denotes a summation.) Substituting Eq. fllOD in 
(PI), we encounter two types of directional terms: nf(PwP)ijfi^ and nfPijfij~. The former 
is nothing but the square of hfPtjVj. One immediately verifies that hf PijuJ = 1 and 
n^PijVj = 0, and therefore 

/- /„s, . 9„ 2-ymB 2 sin 2 9 

h + = 1/2 feo Bin 2 9 = JL ^-z ~ R ~ a ■ 12 

R 1 — p cos 9 

It is illuminating to compare this expression for gravitational perturbations to that of 
the electromagnetic or scalar field of a point source in motion. For a scalar field (f) and 
an electromagnetic four-potential A a in the Lorenz gauge, the standard Lienard-Wiechert 
solution yields (e.g. [18|, p|): 



1 

Q 



— u 7 ■ [x — r (r)]' 



u a (r) 



A a = q — 
r= T0 ' ' -u 7 -[x-r (r)] 7 



(13) 

T=T0 



where q denotes the scalar or electric charge, respectively. Let us denote by the four- 
potential in the transverse gauge. The temporal component Aq vanishes, and the spatial 
part is given by Aj = PijAj. Let A T denote the magnitude of Af , i.e. A T = (AfAf) 1 ^ 2 . 
Then a straightforward calculation yields 

A T = 1 . (14) 

#(l-/3cos0) 1 ; 

To represent all three cases in a single equation, let ip s denote <fi, A T , or h + , for s — 0, 1, 2, 
respectively. Then 
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q (jpsm9) s 

^s = s\-— — , 15 

R 7 (1 — p cost)) 

where in the gravitational case (s = 2) q denotes the particle's mass m. 



III. ANGULAR DISTRIBUTION IN THE ULTRA- RELATIVISTIC CASE 

We shall now consider the angular distribution of the gravitational field in the limit 
7 >> 1, i.e. (3 = 1. We can then omit the factor (3 2 in Eq. ([12]): 

27m sin 2 9 , „ 

h+ = -^^- R n 7 » 1 • 16 

R 1 — p cos 9 

However, the factor f3 in the denominator must be treated more carefully. For 7 >> 1 we 
can always approximate 

1 _ p cos 9 = (1 - 0) + 0(1 - cos 9) = (7~ 2 / 2 ) + (1 - cos 6) . (17) 

Now, for 9 « 1 we can approximate cos 8 = 1 — 9 2 /2, hence 

l-/3cos^= ( 7 - 2 + # 2 )/2 (0«1). (18) 

On the other hand, for # large compared to I/7 we can ignore the first term in the right-hand 
side of Eq. 

l-/3cos#^ l-cos# (0»7~ x ). (19) 

Thus, in the ultra-relativistic limit h + takes two qualitatively different asymptotic forms, 
depending on the value of 9. For small 9 we have 

and for # large compared to I/7, 

/^^(l + cosfl) (^»7" 1 ). (21) 

R 
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The two asymptotic regions overlap at 7 _1 < 5 < 1, where we have h + = A'ym/R. In the 
range of very small angles, 9 < 7 _1 , there is a "hole" in the angular distribution, wherein 
h + sharply decreases and vanishes at 9 = 0. 

Motivated by the beaming phenomenon in the analogous electromagnetic problem, we 
shall refer to the regions 9 < 7 _1 and 9 » 7 _1 as the beaming and off-beaming zones, 
respectively. Note, however, that the metric perturbation h + does not exhibit a true beaming 
phenomenon. There is no sharp enhancement of h + in the narrow forward direction where 
9 ~ 7 ; rather, there is a "hole" at 9 < 7 . (In the off-beaming zone there is a slow increase 
in h + when 9 decreases, like 1 + cos#, but this is a moderate, 7-independent, increase. )[] 

In the ultra-relativistic case, which concerns us here, the beaming zone covers an ex- 
tremely small solid angle, Q oc 7~ 2 . Since there is no enhancement of h + in this range, this 
zone has a negligible contribution to the anticipated rate of detection. Therefore, in what 
follows we shall ignore the "hole" in the beaming zone, and always use the off-beaming ex- 
pression fl2ip for h + . It is remarkable that the perturbation amplitude ([H]) does not depend 
on m or on the particle's speed separately. It only depends on the product 7m, i.e. on the 
particle's energy. 

DtH H analyzed the gravitational field sourced by an ultra-relativistic particle moving 
along a geodesic in flat space. They found that in the limit (3 — > 1 the gravitational field forms 
a planar shock wave which propagates along with the particle, at the speed of light. This 
behavior is very different from what we have found here. This difference cannot be explained 
by the choice of different gauges; For example, in the DtH analysis the geometry before and 
behind the shock wave is strictly flat. The reason for this difference is that, the two analyses 
consider different limiting procedures. DtH took the (constant-speed) boosted Schwarzschild 

1 This qualitative picture changes when one describes the gravitational radiation in terms of the 
Riemann tensor or the effective energy- momentum outflux. In both descriptions, there is a beaming 
effect at angles ~ 7 -1 ; See section V for the energy flux. 
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solution, evaluated the gravitational field at a fixed, finite R, and then took the limit j3 — > 1. 
We are applying here a different limit, which is better adopted to the astrophysical situation 
that motivates the present work: We assume that the particle's speed changes with time, 
then pick the radiative piece of the gravitational field associated with this change of velocity, 
and only then we take the limit (3 — > 1. In this procedure, the non-radiative piece of the 
gravitational field is dropped. (Recall that the TT-projection only respects the radiative 
part of the gravitational field.) This non-radiative piece is unimportant for detection over 
astrophysical distances; yet it is the only piece which exists in the DtH problem, in which 
the gravitational field is non-radiative. 

Let us consider a situation which is perhaps not too realistic from the astrophysical point 
of view, but it may clarify the relation between the two analyses. Assume that a blob is 
ejected from the star at t = 0, with a finite Lorenz factor 7 >> 1, and the observer is located 
at a large R and a very small 9. The blob then moves undisturbed along a geodesic all the 
way from the star to the observer's neighborhood, and passes near the observer. Then we 
expect that, provided that the blob's minimal distance to the observer is small enough (i.e. 
9 is sufficiently small), as the blob passes by, the observer will watch a wave phenomenon 
similar to that described by DtH (though the shock will be somewhat smoothened, due to 
the finite 7). In addition, the observer will also watch the radiative phenomenon associated 
with the ejection, described by the gravitational field fl2l|) . There will be a time lag between 
the two phenomena: The ejection- induced radiation pulse will arrive at the moment t = R, 
corresponding to zero retarded time (we neglect here the small quantity At). On the other 
hand, since the shock-like wave must move along with the particle, it will reach the observer 
at the moment t = R + (i?/2)(7~ 2 - 9 2 ). The two waves will also have different frequencies: 
The typical frequency of the DtH wave will depend on 9 and 7, but will also decrease like 
1/R. The radiative field ([H]) will have an /^-independent typical frequency (which for small 
9 is ~ 7 2 /At; see section V). Therefore, for a detection over astrophysical distances the 
ejection-induced wave seems to be the more important phenomenon. 

The DtH phenomenon will only take place in the range 9 < I/7, because for 9 > I/7 the 
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time lag becomes negative, which would contradict causality. As 9 increases in the range 
< 9 < I/7, the DtH phenomenon will become weaker (because the minimal blob-observer 
distance increases with 9) and the radiative field will become stronger [cf. Eq. (p0|)] . The 
relation between these two gravitational-wave phenomena and their possible coexistence 
deserve further investigation. 

IV. SHOOTING AN ULTRA-RELATIVISTIC BLOB 

We shall now calculate Ah TT , i.e. the change in h (expressed in the TT gauge) that 
occurs when a star ejects an ultra-relativistic blob of matter. As was discussed in section II, 
it is this quantity which is relevant for analyzing the detector's response. We shall show that 
for a sufficiently small blob's energy, this change in the overall metric perturbation coincides 
with the contribution h TT of the blob. 

Consider a star with mass M initially at rest. At a given moment t — to it emits a blob 
of mass m « M with an ultra-relativistic speed 7 >> 1. Let us assume that the whole 
acceleration process ends at t — t\ = to + At (these times all refer to the star's rest frame). 
The mechanism responsible for this process is unimportant for this discussion (it could be, 
for example, an electromagnetic acceleration due to a dynamo effect or some MHD instability 
accelerating a blob of plasma, or radiation pressure). The important point is that at t < t 
the entire system can be modeled as being at rest, then the blob accelerates between t = t 
and t = ti, and at t > t\ the blob moves with a constant speed, 7 >> 1. (One can think of 
this process as the time- reversal of a fully- inelastic collision.) Let us denote the remaining 
mass of the star by M\. We assume that the blob's energy 7m is << M. Then from energy- 
momentum conservation one finds that, at the leading order in the small parameter 7777./M, 
Mi = M — 7m and the star's speed is non-relativistic, (3 S = —■ym/M. 

2 In principle one has to include in this energy-momentum balance the amount of energy E g (and 
also the momentum) carried by the gravitational radiation. However, E g is quadratic in 7m; One 
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Since all motions in this process are along the same axis, the polarization of both the 
blob's and the star's fields will have the same direction. The TT metric perturbation at the 
end of the process can therefore be described in terms of the overall amplitude parameter 
h + (obtained by summing the amplitudes h + of the star and the blob). Correspondingly, 
the change in h TT is given by the quantity Ah + = h l + — h? + , where h? + and h l + denote the 
overall field amplitude before and after the process, respectively. Since initially (3 = 0, from 
Eq. (pD we have h° + = 0. At t > t x we have h\ = hf ar + h b _l° b , where 



h star ^ — ~~~ sin 2 9 , h b i ob -^-(1 + COS 6) . (22) 
R R 

The ratio of these two contributions is 

h star /h Uob MP 2 S sin 2 6 ^ z^ie/Z) 1 — « 1 , (23) 
7m 1 + cos 9 M 

hence the star's contribution to h\ is negligible. We conclude that the change in h + is just 
the contribution of the ultra-relativistic blob: 

Ah + = hl = h b l ob = ^(1 + costf) . (24) 

R 

The fact that Ah + does not vanish means that we are dealing here with a "burst with 
memory" ||. One may be puzzled by this lack of conservation in the overall perturbation 
field, because the source term, the particles' energy-momentum tensor, is conserved. The 
resolution of this puzzle is simple: The Lienard-Wiechert solution (|J) has a dependence on 
the velocity also through the denominator, u 1 [ x — r (r)] 7 . It is this dependence which leads 
to the non- vanishing of Ah + . 

In several astrophysical systems with an accretion disc, the blobs appear to be emitted 
in pairs, along the two polar directions. In such a case, h + will simply be the sum of the 
contributions from the two blobs. Note that the contributions from the two blobs do not 



can show that E g / r ym < 7771/M << 1, so E g may be ignored. The momentum carried by the 
gravitational radiation is also bounded by E g . 
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cancel each other. For example, for a symmetric pair of ultra-relativistic blobs, each carrying 
energy 7m, the sum of the two contributions will be direction-independent: 

Ah + = h\ £ ^ . (25) 

This result agrees with the analysis by Adler and Zeks who considered the case of 
two equal masses [Q. (From the observational point of view, however, recall that the two 
components may arrive the detector at different times, and will also have different Doppler 
factors.) 



V. ANGULAR DEPENDENCE OF OBSERVED BURST FREQUENCY 

Let f c denote the characteristic frequency of the observed burst. It is given by f c ~ 5t~ x , 
where St denotes the burst duration (i.e. the rising time of h + from its initial to final value) 
as measured by the detector. We need to relate 5t to the pulse duration in the star's Lorenz 
frame, At — t\ — t . Let t'(t) denote the arrival time (at the detector) of a null geodesic 
which emerges from the blob at a moment t. These two times are related by t' = t + R. A 
straightforward calculation then yields 

dt' 

^ = l-/3cos0. (26) 
dt 

As was discussed in section III, the beaming zone is extremely narrow and does not signifi- 
cantly contribute to the detection rate. We shall therefore ignore it and use the off-beaming 
approximation (|19|) : 

f^l-cos*, (27) 

and hence 5t = (1 — cos 6) At. Therefore the observed frequency is 

/ c ~ [(l-cos^At]- 1 • (28) 

This relation is independent of 7. Recall, however, that this approximation breaks at 6 < 
7 _1 : At the beaming zone f c saturates at a maximal value /™ ax ~ ^ / At. 
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Although there is no enhancement of h in the beaming zone, the energy flux is in fact 
beamed in the forward direction, as we now show. Let F denote the time-integrated energy 
flux per unit solid angle. Since the energy density is proportional to h 2 uj 2 ~ h 2 /St 2 , we have 
F ~ R 2 h 2 /5t ~ ( , ym) 2 /5t. Thus, outside the beaming zone (but at 9 « 1) F behaves as 
[(jm) 2 At -1 ] 6~ 2 , and it gets a maximal value of order / y 2 ( , ym) 2 At~ 1 at 9 ~ I/7. 

This result is remarkable, because no such beaming occurs in the analogous situation of 
an ultra-relativistic head-on collision of two black holes [ 14[] , [16| . To sharpen the contrast 
between the two cases, consider the time-reversal variant of our problem, i.e. a blob colliding 
fully inelastically with a star. As one can easily verify, the pattern of Ah in this collision 
problem will be exactly the same as in the original ejection problem - and, in particular, 
F will be beamed in the forward direction. The intuitive reason for the difference between 
the two collision problems is simple: In our star-blob system, the gravitational field is 
everywhere weak; the interaction between the star and the blob is non-gravitational. In 
particular, the (de-) acceleration occurs on a distance scale At which is >> M. Hence the 
beamed gravitational energy flux propagates to null infinity without any obstacle. On the 
other hand, in the analogous black-holes head-on collision problem the interaction is solely 
gravitational, and the (de-) acceleration occurs on a distance scale which (to the extent it is 
defined) is of order M . Hence, should any beamed radiation form, it would immediately be 
swallowed by the large black hole. It is therefore not surprising that no beamed radiation is 
observed at null infinity. [^1|] Indeed, in the case of an ultra-relativistic black-holes encounter 



with a large impact parameter, the gravitational radiation does exhibit a beaming [10], [14 



3 This argument does not imply that a beamed radiation will actually hit the large black hole: 
Since the deacceleration distance in this case is of the same order of magnitude as the radius of 
curvature of the large black hole, M, the qualitative flat-space arguments are not applicable for 
the analysis of the behavior of radiation near the blach hole. In fact, the equivalence principle 
suggests that no beamed radiation will hit the black hole. 
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15], [16 . 



VI. CALCULATION OF DETECTION RATE 



In this section we shall evaluate the burst's detection rate, as a function of the various 
parameters involved. To this end, we shall first consider the detector's sensitivity. Then 
we calculate, for a given burst, the detection distance, the detection angular range, and the 
detection volume. Given the rate of such events of blob ejection, we shall derive the general 
expression for the detection rate. 



Let us denote the detector's peak sensitivity by h d} and the frequency at which this 
maximal sensitivity is achieved by fa. For a burst with memory, the detector will have the 
maximal sensitivity h d as long as the burst's observed frequency f c is larger than f d ||. The 
detector's sensitivity is quickly degraded at fa < fa, and in the calculation below we shall 
neglect this range for simplicity (this may result in a small decrease in the calculated detec- 
tion rate). Thus, we shall presume that a detection occurs if the following two conditions 
are satisfied: 



Detector sensitivity 



h + > h d 



(29) 



and 



fa> fa ■ 



(30) 



Detection distance 



From condition (|29|) and Eq. (|2l|), at a given direction the observation distance is 




(31) 
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The maximal observation distance is obtained at the forward direction, 

#max = 47m/ h d . 



(32) 



Detection angle 



From condition (|30|) and Eq. (|28|) , the angular detection range is bounded by 



(1-COS0) < (f d At) 



(33) 



We shall distinguish between two cases: 

Case A- At < 1/(2 fa): In this case, the burst may be observed in all directions. 

Case B- At > l/(2f d ): In this case, the burst is observed in the range 6 < max , where 
#max is given by 



There is a third case, in which At is of order 7 2 //d or larger. In this range, the observed 
frequency will be too small even at the forward direction. Hereafter we shall assume At << 
l 2 / fd ■ (This assumption seems very reasonable if we take e.g. 7 ~ 10 3 ; See the discussion 
of the values of the various parameters in the next section.) Note also that for such a large 
At, # max of case B is so small that the resulting detection rate is negligible anyway. 



Consider an event of a blob ejection. We shall now calculate the volume of the region of 
space in which the emitted burst will be detected. This volume is given by 



1 — cos 6 



max 



(f d At) 



(34) 



Detection volume 



1 




(35) 
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(In case A, we take 6* max = tt.) Substitution of Eq. ( [31]) in the last expression yields 

V = ^f( 7 m/h d ) 3 w , (36) 

where 

w = l-(l + cos^ max ) 4 /16 . (37) 

Evaluating w in the two cases A,B, we find 

f 1 (At < At d ) , 

io = I (38) 

[ 1 - (1 - At,/At) 4 (At > At d ) , 
where At d = 1/(2 fa). In the case of large At, we get 

w = AAt d /At (At » At d ) . (39) 

Detection rate 

Let us denote by n the rate of supernovae explosions per unit volume. We assume here 
that each supernova ejects blobs to the two polar directions, which adds another factor of 
two {] The detection rate is therefore 

N = 2nV =^^n(^m/h d f w . (40) 
3 

According to the cannonball model, in each of the two polar directions several blobs are 
emitted, typically of the order 3-5. We do not multiply 2nV by this number, because these 
are not independent detections. Rather, each observed event will be a composition of a few 
bursts.PI 



4 This factor 2 is not mathematically precise, and it may depend on At/ Atd, but in the worst 
case - case A - it is 15/8, so we can well approximate it by 2. 

5 It is assumed here that the time separation between two successive bursts in a given event will be 
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VII. INSERTING ASTROPHYSICAL PARAMETERS 



We shall now evaluate the detection rate by substituting astrophysical numbers in the 
general expression, Eq. fl4*0|). The astrophysical situation concerned us here is that of a 
supernova explosion resulting in the ejection of ultra-relativistic blobs, as proposed by the 
cannonball model. First we re-write the above general expression, retaining the constants 
C and G: 



Ar 128tt fymGY / x \ 



where w is given in Eq. ([38]). We need to evaluate the supernova rate n, the blob's energy 



7m, and the characteristic acceleration time At. 

Consider first the supernova rate n. The Shapley-Ames 'fiducial' sample of 342 galaxies 
within the Virgo circle |[22|| , [p3fl has a mean B-band luminosity of 6.7/i~ 2 10 9 Le (B) and a 
supernova rate of 3.09/i 2 [lOO?/rlO lo L0 (B)] 1 . The luminosity density of the local universe 
g4[ p L = (2.0 ± 0.4) • 10 8 hL e Mpc- 3 . Therefore we estimate the supernova density in the 



is 



local universe (taking h = 0.65) as 

n ~ 1.7 ■ lO^Vr^ 1 ■ Mpc' 3 . (42) 

For the blob's energy, the cannonball model suggests the typical value E = ^mc 2 ~ 
10 52 erg 0, corresponding to 7m of order 10 31 gram. Using this and Eq. (|42|) one gets 

"7.5- 10- 23 ! 3 



N 



where E 52 is the blob's energy in units of 10 52 erg. 



El 2 w(At) [Yr~ l ] , (43) 



large compared to Atd, so they will not sum up coherently (a coherent sum-up would significantly 
increase the detection distance and hence the detection rate). Note, however, that the presence of 
several bursts in a single event eases its detection: It makes it easier to distinguish the event from 
the noise, so it may decrease the signal-to-noise ratio required for detection. 
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The sensitivity curve for LIGO-II (advanced detector) may be found in e.g. Ref. [25 . 
The optimal value of the noise level is about 6 ■ 10~ 24 , at a frequency f d ~ 50 sec" 1 . A 
minimal signal-to-noise ratio of about 5.5 |26| is required for a detection. [|27| An averaging 
over all possible source polarizations yields an extra factor ~ 0.5 in the effective wave's 
amplitude. Combining these two factors, one obtains an effective threshold for detection 
which is about 11 times larger than the noise level H|, i.e. h d « 7 • 10~ 23 . This leads to 



iV ~ 1.3 E£ 2 w(At) [Yr- 1 ] . (44) 

The characteristic time At is hard to estimate, because the mechanism responsible for the 
acceleration is unknown (presumably it is the same yet-unclear mechanism which forms the 
relativistic jets in various astrophysical systems including accretion discs around a compact 
object). The parameter w depends on the dimensionless parameter At/At,, where for 
LIGO-II Atd ~ 10~ 2 sec. As long as At is smaller than At, we have w = 1 and hence 

iV~ 1.3 E^Yr- 1 ] (At<At d ) . (45) 

For At > Atd, w decreases. If At >> At d we may use the approximation ( p9|) and the 
detection rate is inversely proportional to At: 

N ~ 5(At d /At) E^Yr- 1 ] (At » At,) . (46) 

A reasonable lower bound on At may be obtained by assuming that the acceleration to 
an ultra-relativistic speed occurs within a distance scale comparable to the neutron star's 
radius, R ns ~ 10km, namely At min « R ns /c ~ 3 • 10~ 5 sec. 

Dar p8| suggested that the typical time scale for the whole blob ejection process should be 



comparable to the dynamical gravitational time-scale, e.g. the free-fall time at the neutron 
star's radius. This time scale is of order Atf ree _f a u ~ 10~ 4 sec. According to this suggestion 
we can safely use the detection rate (|5|). Since Atf ree _f a u « At d , this conclusion will 
hold even if At is larger than Atj ree _j a ;/ by two orders of magnitude. 

To conclude, according to the cannonball model, a reasonable estimate of detection rate 
by LIGO-II is about Ef j2 per year (assuming a small At). The parameter E 52 is estimated 
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to be ~ 1, but there is an uncertainty of almost an order of magnitude. This leads to an 
uncertainty of order ~ 10 2 in the detection rate. The latter will be smaller if At is larger 
than Atd- 

The maximal observation distance, achieved in the forward direction, is given in Eq. 
(|52]). Substituting the above astrophysical parameters, we obtain for the advanced detector 

#max = 4^ « 15 • E b2 Mpc . (47) 
c 2 h d 

The cosmological redshift is negligible at this distance. 
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